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Abstract 

A nonlocal energy-balance equation is derived for the optical absorption, photoluminescence and 
inelastic electron-phonon scattering, which determines the electron and hole temperatures for any 
given lattice temperature. The evolution of the lattice temperature is found to be determined by 
the difference between the power-loss density due to photoluminescence and the power-gain density 
due to optical absorption, as well as by the initial lattice temperature. We find that in addition 
to the expected decrease in the lattice temperature, the electron temperature also decreases with 
time. A laser-cooling power as high as 380 eV/s is predicted for the wide bandgap semiconductor 
A1N initially at room temperature when the pump-laser field is only 10 V/cm. Laser cooling is 
found to be more efficient for a large bandgap material, a weaker laser field, and a high initial lattice 
temperature. The laser-cooling rate then decreases as the lattice cools. The theory presented here 
provides quantitative predictions that can guide future experiments. 
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The coolingof a solid via light-induced fluoresce has been of interest for a very loner 
This interesting phenomenon involves the excitation of an electron 
from the valence bandedge to the conduction bandedge by absorbing a pump photon. This 
cool electron quickly becomes hot by gaining thermal energy through ultrafast electron- 
phonon scattering. After a radiative lifetime, recombination of the hot electron will produce 
a spontaneous photon with energy higher than that of the pump photon. As a result, the 
lattice will be cooled due to the loss of thermal energy to the electron. It is only recently 
that this phenomenon has been observed experimentally. Laser-induced fluorescent cooling 
of heavy-metal-fluoride glass doped with trivalent ytterbium ions was the first realization 
of this concept jfj. Soon to follow were demonstrations of cooling in dye solutions |?| and 
thulium-doped glass 

Laser cooling of a semiconductor, however, remains an elusive goal, although it has been 



pursued for many years 



ft 



101 ] . Indeed, now more than ever the field of laser cooling 



is a topic of intense on-going theoretical and experimental investigations 11]. The only 
theories attempting to model the laser cooling phenomenon are local simulation theories 
that include rate equations for determining the steady-state carrier density and the loss 
of lattice energy with several kinetic coefficients. The problem with these theories is that 
they neglect important dynamical effects such as the change of the carrier distribution when 
the temperature is lowered. Therefore, they only apply to situations with little change 
of temperature. The main feature of the rate equation approach is its simplicity, but it 
is unable to elucidate the essential physics behind the laser cooling phenomena. The key 
question that still remains open is, what are the best semiconductor materials and conditions 
for achieving the greatest laser cooling effect? This requires an accurate nonlocal theory on a 
microscopic level, which directly provides an evolution equation for the lattice temperature 
by including the dynamical effects. This theory should establish a criterion for the occurrence 
of laser cooling. The theory should also establish a criterion for the efficiency of the laser 
cooling if it does exist. In this paper, we will present such a nonlocal theory for the laser 
cooling of semiconductors. By including the effect of the carrier distribution, we will be 
able to uncover the essential physics underlying the phenomenon and we will be able to 
provide important quantitative predictions that can guide experimentalists toward achieving 
maximum efficiency of laser cooling in the future. 
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Let us first recall the transport force-balance theory for hot electron transport Q|. When 
a dc electric field is applied to a doped semiconductor, there is a drift of electrons as a 
result of the center-of-mass motion. This is described by a balancing of forces between the 
frictional forces due to scattering processes and the electrostatic force. At the same time, 
electrons form a quasi-equilibrium Fermi-Dirac distribution. The electron temperature of 
this distribution can be determined by the energy balance between the power-gain density 
due to joule heating and the power-loss density due to electron-phonon scattering. As a 
result, the electron temperature becomes higher than the lattice temperature if the lattice 
is in thermal equilibrium with a heat bath and the drift velocity is large. 

For the situation under consideration in this paper, a weak pump laser first excites elec- 
trons from the valence bandedge to the conduction bandedge. The excited carriers instanta- 
neously form a non-equilibrium distribution [3]. It is well known that the quantum kinetics 
of the scattering of electrons with phonons or other carriers under a weak pump field can 
only be seen within the time scale of several hundred femtoseconds [14]. Subsequently, ul- 
trafast carrier-phonon and carrier-carrier scattering quickly adjusts the kinetic energies of 
these excited carriers by taking energy from the lattice Q|. As a result, a quasi-equilibrium 
Fermi-Dirac distribution of carriers is 

ture determined by the pump-field intensity, pump-photon energy, and lattice temperature. 
After a few tens of nanoseconds, radiative decay of the excited carriers will begin to af- 
fect the electron distribution. The electron temperature will be adiabatically readjusted 
according to an energy balance between the power-gain density due to optical absorption, 
the power-loss density due to photoluminescence, and the power-exchange density due to 
scattering with phonons. At the same time, the lattice temperature will evolve because 
of an imbalance between the power loss due to transferring phonon energy to carriers and 
the slight power gain from the external thermal radiation. Just before the radiative decay 
occurs, the lattice and the electrons are in thermal equilibrium with an initial temperature 
which can be determined by solving a semiconductor Bloch equation [lfi| . 

Using the nonlocal theory described below, we find that the laser-cooling rate is largest for 
a large bandgap material, a weaker pump-laser field, and a high initial lattice temperature. 
We also find that the laser-cooling power decreases as the lattice cools down. 

In general, the power gain by electrons from the absorption of the pump laser can not 
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be balanced by the power loss due to spontaneous photon emission alone. As a result, 
electrons either take energy from or give energy to phonons through inelastic scattering, 
which depends on the sign of the difference between the electron and lattice temperatures. 
The electron temperature can be determined by an energy-balance equation for any given 
lattice temperature. The pair scattering between electrons due to the Coulomb interaction 
conserves the total energy and does not contribute to the energy-balance equation. On the 
other hand, the single-particle electron-phonon scattering greatly contributes to the energy- 
balance equation. 



The power-density loss due to spontaneous photon emission is calculated to be [16, 11 
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where /| (/£) is the quasi-equilibrium electron (hole) distribution at electron (hole) temper- 
ature T e (T h ), E G is the bandgap, A is the spin-orbit splitting, m* is the electron effective 
mass, fi* is the reduced mass of electrons and holes, e r is the average relative dielectric con- 
stant, V is the volume of the sample, and E' G (k) is the renormalized bandgap. It is clear 
from Eq. Q that the larger the bandgap or the higher the carrier temperature, the stronger 
the power-loss density will be. 

The power-density gain due to pumping by a spatially-uniform laser for hQ p > E G is 
found to be [l8| 
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(2) 



where 70 is the homogeneous level broadening, Ml p is the pump-photon energy, and 2Afe 
is the renormalized Rabi splitting which is proportional to the square root of the pump- 
laser intensity. It is seen from Eq. (j2J) that the greater the pump-laser field, the higher the 
power-gain density will be. By comparing Eqs. Q and (0), we find that the ratio of the 
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power- loss density to the power-gain density scales as Eq/S^. For simplicity, we do not 
include recapture of photoluminescence photons here, which can be equivalently included as 
an adjustment of the pump-laser intensity if the carrier temperature is much smaller than 
E G /k B . 

By keeping only the leading order interaction between electrons and phonons or impurities 
in the Heisenberg equation, we get the following power-exchange density W s e from impurities, 



phonons and scattering-assisted photons to electrons 
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which includes the phonon- or impurity-assisted photon absorption for n ^ 0. Here, £ p is 
the pump-laser strength, J5| is the renormalized electron kinetic energy, and T e and Tl are 
the electron and lattice temperatures, q\\ lies in the polarization direction of the pump-laser 
field, N^ h (x) = [exp(x) — is the Bose-Einstein function, n is an integer, J n (x) is the nth 
order Bessel function, hu q x is the phonon energy for phonon wave number q and mode A, 
rii is the impurity concentration, \U?(q)\ 2 = |e 2 /[eoe r (g 2 + A 2 )V]| 2 is the electron-impurity 
coupling strength, 1/A e is the static screening length, and | A | 2 is the electron-phonon 
coupling strength. For polar semiconductors, such as A^Gai^N, there exist both acoustic 
and optical phonon modes. For optical phonon modes, only the longitudinal-optical phonon 
mode can strongly couple to electrons. For acoustic phonon scattering, on the other hand, 
we use the deformation-potential approximation 2l| with parameters given in the text. The 
detailed form of electron-phonon coupling strength |C| A | 2 can been found from publications 
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22l | . Applying the Debye model to low-energy acoustic phonons, we get u q \ = c\q 
with X = £, t. For holes we get the similar power-exchange density W s h with It is 
clear from Eq. (JHJ that the electron energy loss or gain from phonons under weak pump field 
depends on whether the electron temperature is higher or lower than the lattice temperature, 
respectively. The same argument applies to holes. 



In order to cool the lattice, the power gain of the electrons due to laser pumping must 
be smaller than the power loss due to spontaneous photon emission. This requires a very 
weak laser field and a large bandgap. The energy conservation in steady state requires 
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W ab - W sp + W s c + W s h = . (4) 

The solution of this equation provides the carrier temperature for any given lattice tem- 
perature T L since Wg + explicitly depend on the lattice temperature T L . The sign of 
W a b — VV sp determines the signs of T e — Tl and 21 — Tl. The larger the value of | W a b — VV sp | , 
the larger the deviation of carrier temperature from Tl will be. 

Although the phonons also stay in a quasi-equilibrium state, the phonon temperature 
Tl directly evolves with time due to an imbalance between the power loss to electrons and 
holes and the power gain from any external thermal source (such as the background thermal 
radiation). As a result, the average phonon energy varies with time. This gives rise to 

v (T B 4 - Tl) - (W sp - W ab ) , (5) 

where ol is the lattice constant, a = 7i 2 k^/60h 3 c 2 is the Stefan-Boltzmann constant, A, is the 
surface area of the sample. We assume T L = T at t = 0, where T is the initial temperature 
of equilibrium phonons, and Tb is the environmental temperature which is close to To for 
bandedge pumping with very weak laser field. The first term in Eq. © is much smaller than 
the second term even when Tl 7^ Tb. The rate of reduction of Tl is determined by Wg + Wg 
which decreases with decreasing Tl and the temperature difference Tl — T c . Moreover, we 
note that the bandgap of semiconductors in general depends on the lattice temperature, 
but can be neglected for wide bandgap semiconductors such as Al^Gai.^N. From Eq. (jSJ) 
we know that the cooling of the lattice implies W| + W s h = W sp — W a b > 0. This requires 
T c < Tl (Th < Tl for hole) from Eq. (jHJ) with a weak pump field and a large bandgap. 

In this paper, we consider the semiconductor Al x Ga\- x N for our numerical calculations, 
where x is the percentage of Al in the alloy. The bandgap increases with x. 

For Al x Gai_ x N ', we choose the following parameters: Eq = 3.4 + 2.8x eV, m*/mo = 
0.2 + 0.2s, ml/m = 1.4 + 2.13x, A = 0.02 - O.OOlx eV, e s = 8.9 - 0.4a:, = 5.35 - 0.58a;, 
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FIG. 1: Calculated scaled lattice temperature T^/Tq for Al^Gai-^N as a function of time t with 
M7 p — Eq = 10 meV for four different cases. These cases include: x = 1 and £ p = 10 V/cm, 
To = 300 K (solid curve); x = 0.5 and £ p = 10 V/cm, To = 300 K (dashed curve); x = 1 and 
£ p = 100 V/cm, T Q = 300 K (dash-dotted curve); and x = 1, S p = 10 V/cm, T = 40 K (dotted 
curve). The other parameters are given in the text. 

hu LO = 91.2 + 8.0x meV, p = 6.15 - 2.92x g/cm 3 , D = -(8.3 + 1.2x) eV, h u = (2.81 + 
4.09a;) x 10 7 V/cm, c t = (2.68 + 1.02a;) x 10 5 cm/sec, c e = (6.56 + 2.56a;) x 10 5 cm/sec, 
a L = 5.12 - 0.14a;, 7 = h/r with r = 0.1 ps, n { = 10 10 cm" 3 , e r = (e s + e 00 )/2, hQ p - E G = 
10 meV, and the sample is assumed to be cubic with an edge size of 1 cm. 

Figured displays our main results for the scaled lattice temperature Tl/T as a function 
of time t (0 < t < 1 /as) for Al^Gai-^N. From it we find that the laser cooling at To = 300 K, 
x = 1.0 and S p = 10 V/cm (solid curve) is the largest compared to the other three cases, 
reaching as high as k B AT L /At = k B (T — T L )/ At = 380 eV/s. The cooling effect becomes 
smaller when x is reduced to 0.5 (dashed curve) with a smaller bandgap. Moreover, the 
laser cooling changes into laser heating when S p is increased to 100 V/cm (dash-dotted 
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curve). Finally, the laser-cooling efficiency decreases to 100 eV/s when T drops to 40 K 
(dotted curve). This indicates that the laser cooling of a lattice can be maximized for wide- 
bandgap semiconductors under the conditions of low pump-field strength and high initial 
lattice temperature. 

In conclusion, by using the energy-balance equation for pump-laser induced conduction 
electrons and holes, we have demonstrated a laser-cooling power as high as 380 eV/s for 
the wide bandgap semiconductor A1N at room temperature when the pump-laser field is 
only 10 V/cm. The evolution of the lattice temperature was found to be determined by 
the difference between the power-loss density due to photoluminescence and the power-gain 
density due to optical absorption, as well as the initial lattice temperature. 

The authors are grateful for helpful discussions with Prof. M. Sheik-Bahae from the 
University of New Mexico, Dr. B. Flake from the Air Force Research Lab, and Dr. R. I. 
Epstein from the Los Alamos National Laboratories. 
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